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Introduction 

Body-conforming coordinates transformation of a fluid conservation-law form is 
generally used in computational fluid dynamics. The metrics associated with the coordi- 
nates transformation are required to satisfy certain geometric identities to maintain the free 
stream. 1 These metrics are called free-stream capturing (or preserving) metrics. So far, nu- 
merical techniques are known to capture the free-stream on stationary grids. 2-4 However, 
the extension of the free-stream capturing metrics to moving grids is not straightforward. 
The error introduced by the time metrics has been overlooked because it is negligible in 
most cases, but it can be significant in certain applications such as helicopter rotor flow 
fields. 5 

Rigorous formulations to avoid this error were su gg ested in Ref. 1, and demon- 
strated, for example, in Ref. 6. Based on the work in Ref. 1, the present study describes 
detailed formulas for constructing the free-stream capturing metrics in space and time on 
both the finite-volume (FV) and finite-difference (FD) framework. The error introduced 
by the inconsistent time-metric term is also evaluated. 


^Research Scientist, Applied Computational Fluids Branch, MCAT Institute, San 
Jose, California 95127. 
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Finite- Volume Formulation 


Geometric Identities and Free-Stream Capturing 

Following Ref. 1, the integral form of a conservation law for a given cell can be 
written as 

[ QdV — f QdV+ [' <f n • FdS dt = 0 (1) 

Jv(ti) JV(t x ) Jti JS(t) 

where V(i) is the cell volume and n dS(t) is a vector element of surface area with outward 
normal n. Considering the Euler equations, Q is a vector of conservative variables, viz., 
density, momentum and energy, and F is the flux tensor of Q. The flux F can be de- 
composed into the flux in the stationary frame, F J( , and the contribution due to surface 
element velocity, v as follows 

F = F 3t - vQ (2) 


Let r, r 0 (t), v 0 (t), and il(t) be the position vectors of a point in space, the origin, velocity 
and angular velocity of the non-inertial frame relative to the inertial frame, respectively. 
Then, 


V = V r + V c 


(3) 


where 


Vr = V 0 (t) + X [r - r 0 (t)] 


(4) 


and v c (<) is the surface element velocity relative to the non-inertial frame, but expressed 
in the inertial frame. 

The geometric identities given in Ref. 1 are as follows. A mathematical expression 
of a closed cell is given by 

j ndS = 0 (5) 

The relative rigid motion of two frames of reference is given by 


£ n • v r dS = 0 


( 6 ) 
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Using Eqs. (4) and (6), one obtains 


£ r x ndS = 0 


(7) 


The conservation of volume for a time- varying cell is given by 


V(t2)-V(ti)= f <£ n -\ c dS dt (8) 

Jh Js(t) 

The free-stream preservation due to these geometric identities can be demonstrated by 
substituting Qoo, F s too Eqs. (2) to (4) to Eq. (1): 

Qoo[V(t 2 ) - V(<i)] = - £ \F 3too ■ ( jf n dS) - QooCvo - Q x r 0 ) • (jf ndS ) 


- Qoott ■ r X n dS) - Qoo(j> n • v c dS)]dt 


(9) 


The geometric identity, Eq. (5), suffices to capture the free stream in a fixed coordinate 
system, where most of steady-state computations are carried out, and in a moving coor- 
dinate system without rotation (Q = 0). When the grid is moving with rotation (fi ^ 0), 
the second geometric identity, Eq. (7), is to be satisfied. For the general motion of grid 
with changing cell volume, the third geometric identity, Eq. (8), is also required. 

The geometric identity, Eq. (5), preserves the free stream at any instance t when 
v = 0. Thus, the time-differential form of Eq. (1) is often used in the FV formulation. 
However, if the grid moves, the geometric identities have to be satisfied correctly in the 
integral form. 


Free-Stream Capturing in the Inertial Frame 

To preserve the free-stream perfectly with a moving grid, Ref. 1 suggests to consider 
the rigorous FV formulations in space and time. One of the rigorous FV formulations with 
a grid velocity expressed in the inertial frame is shown here in detail. 


Assuming that v is given and then by substituting Q <*> and T 3 too to Eq. (1), one 


obtains 


Qoo[V(t 2 ) - V{ti)] = -F 3too . jf ^ndSdt + Qoo jfn -wdSdt (10) 
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The first term in the right-hand side of the above equation yields Eq. (5). Then the rest 
of the equation becomes similar to Eq. (8). 

Figure 1 shows a regular hexahedral cell. We assume that all edges are straight 
lines. Reference 1 describes the formulas for the surface vector S (note that S is taken 
in the positive coordinate direction here) so as to satisfy the geometric identities on the 
hexahedron: 


Si562 = ^(r 6 - ri) x (r 5 - r 2 ) 

(11) 

Si562 = (r56 - ri 2 ) X (ri5 - T 26 ) 

(12) 

1^12345678 = g(Sl485 + Sl234 + Si56 2 ) * (f7 _ ) 

(13) 


where r 56 = |(r 5 + r 6 ), and so on. Note that Eqs. (11) and (12) result in the same 
expression. In fact, the surface vector is defined uniquely as long as the edges are straight 
lines. Also note that there are other consistent ways to compute the cell volume instead 
of Eq. (13). However, Eq. (13) is the simplest form. Let S 1562 == >^i562n. With either 
Eq. (11) or (12), one obtains 

]£Sn = 0 (14) 

cell 


Thus, Eq. (5) is satisfied. 

It is essential to compute the second term in the right-hand side of Eq. (10). It can 
be rewritten as 

f (^ Sn • v)dt = ^( / Sn-vdf) (15) 

•'*1 cell cell *' <l 

Let Vs be a volume swept by a surface S between the time interval [<i , £ 2 ]- 


Vs 


= /" 

*/<! 


Sn • \dt 


(16) 


Let S(tj) = 5 is 62 and S(t 2 ) = Svs>e>2> (see Fig. 2). The volume V s can be computed 
similarly to Eq. (13) as 

VS 1234 =V 7 122 ' 1 ' 566 ' 5 ' 


= -(Sll' 5'5 + Si 22 'l' + S 1562 ) • (r 6 < — Fl ) 


(17) 
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where Sn'5'5 an d Si 22 'i' are the surface vectors in space and time domain. Note that this 
formula requires only the difference of the positions of grid points between 1 1 and <2> not 
the grid velocity v itself. Then, instead of Eqs. ( 6 ) and ( 8 ), one obtains 

AV = Tv s ( 18 ) 

cell 

This identity does not mean to satisfy Eq. ( 7 ) but does satisfy Eq. ( 6 ) in the time- 
integral form. Therefore, Eq. ( 18 ) leads to the perfect free-stream capturing with the 
use of Eq. ( 11 ). 

Free-Stream Capturing in the Non-Inertial Frame 

It is convenient to use the non-inertial frame for certain applications. Thus, the 
FV formulation with a grid velocity expressed in the non-inertial frame is shown next. 
The analysis also provides a deep insight for the free-stream capturing because it considers 
three types of motion given in Eqs. ( 3 ) and ( 4 ) separately. The discretized forms of the 
geometric identities in the FV method can be expressed as 


M 

to 

3 

II 

0 

(19) 

cell 

VS r x n = 0 

(20) 

cell 

AV = j2v Sc 

(21) 


cell 


where Vs c is obtained from Eq. ( 16 ) by replacing v with v c (see Ref. 1 for more details). 

Reference 1 introduces the area moment M = J s r x ndS so as to satisfy the 
discretized geometric identity, Eq. ( 7 ), on the hexahedron: 

Ml562 = Tl65 X Sl65 + 1*126 X Sl26 (22) 

where nes = 3(1*1 + r 6 + r 5 ), Si 6 s = 3(1*6 - iq) x (r 5 - ri), and so on. Note that 
M1562 7^ ri562 x Si562- The expression, ri562 x S1562, is not well-defined for computing 
area moment. In contrast, Eq. ( 22 ) is well-defined. To see these, let r^, r^, S<> and Sa 
be ri562, ri65, |Si562|n and |Si 6 s|n, respectively. After simple algebraic manipulation, one 


5 



obtains 


r 0 X S<> 7^ 0 


(23) 


cell 


in contrast, 


53 r A X S A = 0 


(24) 


cell 


Therefore, in addition to Eq. (21), the free-stream capturing can be shown as 

[F 3 <oo - (v 0 - Cl x r 0 )<2oo] • (53 S <>) “ ‘ (53 M ) = 0 ( 25 ) 

cell cell 

The inconsistency of Eqs. (23) and (24) was pointed out in Ref. 1 and had been overlooked 
for constructing the free-stream capturing metrics in space and time. In other words, 
surface vectors and moments have to be computed to satisfy Eqs. (5) and (7), respectively. 

The error introduced by the use of the inconsistent area moment has a unique 
feature. Let the cell surface be a parallelogram (S<> = S A + S v = 2S A and r<> = 
i(r A + *v)). Now the difference between Eqs. (23) and (24) can be shown as 

53 (*■$• x S A - r A x S a - r v x S = 53 ([( r 0 ” r A) + ( r 0 ~ r v)] x = 0 (26) 

cell ce M 

Thus the error disappears, for example, on the Cartesian grid. For a hexahedron having 
arbitrary quadrilateral surfaces, the error can be written as 

^l53(r 0 -r A )xS A | = 0( 1 ) (27) 

cell 

where |r<> -r A | = 0(h), |S A | = 0(h 2 ) and 7 = 0{h 3 ). This error may be ignored as long 
as the effect of Coriolis force is negligible in the flow field. 


Alternate Ways of Free-Strearn Capturing 

There Eire two other ideas to preserve the free stream. The simple way, especially 
for the rigid motion of a grid, is to use the free-stream subtraction technique, 2 when the 
free stream is uniform. One can obtain the equation by replacing Q and F in Eq. (1) with 
Q - Q 0 o and F - Foo, respectively. 
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The other conceptually different way is to use Sa always; that is, to regard the 
hexahedral cell as dodecahedron or to divide the hexahedral cell into tetrahedron. Then, 
instead of Eq. (14), one will obtain 

(28) 

cell 

for either dodecahedron or tetrahedron in addition to Eq. (24). The resulting metric terms 
will preserve the free stream. The use of the tetrahedral cell allows the most compact and 
consistent metric formulation. Note, however, that the use of the tetrahedron results in 
unstructured-grid formulations. 

Finite-Difference Formulation 
Geometric Identities in the Finite-Difference Formulation 

The analysis of the FD method can be simplified with the aid of the above discussion 
of the FV method. The FD formulation has to be derived from the integral form, Eq. (1). 
Again from Ref. 1, the differential form for Eq. (1) can be written with a generalized 
coordinate transformation, 

r = r(£,7?,C,T), t — t (29) 

as follows: 

Q r + + F„ + G< = 0 (30) 

where subscripts indicate partial differentiation, 

Q = QV, E = S* • F, F = S" • F, G = S c - F 

and where 

S ? = r,, x r<;, S’* = r<; x r^, S'* = r { x r^ (31) 

and 

^ = r{-r,xr ( ; (32) 
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These metrics axe related with the usual FD notations as 


and 


Yi 

j 




(33) 



(34) 


where r r = v and J is the transformation Jacobian, J — 1/V. Analogous definitions can 
be derived for the other directions. The differential forms of the geometric identities are 
known as 

(S«)« + (S ,? ), + (S <: ) c = 0 (35) 


and 

V r = (S f ■ r,)< + (S’ ■ r,)„ + (S< ■ r r )< (36) 

Equation (36) is called the differential statement of the geometric conservation law (GCL). 7 

FVee-Stream Capturing Metrics in Space and Time 

Following Ref. 1, let the edges of the hexahedron in Fig. 1 be redefined as a double- 
sized cell in the FD grid (ri = r»— i,j— i 1 , r 2 = ri-ij+i.t-i, r 3 = r.-ij+i^+i, • • •, 
r 8 = n+ij-i'k+i)- Also let the time level advance from t\ to < 2 - Then, all the discussions 
for the FV formulation in the previous section can be applied to the FD formulation for 
the central differencing. 

The surface vector evaluations, Eqs. (11) and (12), can be regarded as the evalua- 
tions of the free-stream capturing metrics in the stationary grid for the FD method. For 
example, Eq. (11) has been applied for the central-differencing part in Ref. 4. Equation (12) 

can be rewritten as 

s Li,i,k = ( r 56 - ri 2 ) x (ris - r 26 ) 

(rs - ri) + (r 6 - r 2 ) w (ri - r 2 ) + (r 5 - r 6 ) 

= 2 2 

— — ~(^c r i-i,j-i,* + <5<Ti_ij+i,t) x + 6 v r,_ij,i+i) 
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where 8 denotes the difference operator. After scaling by one-fourth to adjust the area from 
double-sized to regular cell, the first component of the above expression can be written as 

y = (^y)(Mc 2 ) - VXMcv) ( 3S ) 

where n denotes the arithmetic averaging operator. Thus, Eq. (12) is equivalent to the 
consistently differenced metrics in Ref. 3 that are based on the averaging procedure so as 
to satisfy the differential chain rules numerically. 

The main discrepancy between the FV and FD formulations appears in the defini- 
tion of cell volume. The cell volume defined by Eq. (13) is different from the one defined 
by the discretized form of Eq. (32) because the FD method does not use the cell concept. 
Nevertheless, Eq. (13) can be applied to the FD method with a scaling factor of one-eighth 
instead of Eq. (32). Then, the FV space metrics on the double-sized cell become identical 
to the FD ones. 

The FD time-metric evaluation is also considered from the FV point of view. It is 
easily found that the time-metric evaluation, Eq. (34), will not maintain the free stream 
even with the use of the free-stream capturing metrics, Eqs. (11) or (12), in case of a 
rotating frame because of Eq. (23). Also, it can be shown that such inconsistent time 
metrics do not satisfy GCL. The discretized form of GCL can be written as 

AV = Arfo(S< • r r ) + ^(S" • r r ) + 6 C (S< ■ r r )] (39) 

Let the grid move in the rigid rotation, that is, V T = 0 and rr = fixr. Then the left-hand 
side of Eq. (39) is zero. But the right-hand side results in ^(rxS^)+^(rxS' , )-|-^(rxS^) ^ 
0 (Eq. (23) appears again). This indicates that the use of the GCL condition, Eq. (39), 
for computing AV can be erroneous. In other words, the GCL condition, Eq. (39), is 
necessary to preserve the free stream, but not sufficient to construct consistent metrics in 
space and time. 

It is easily found that the time integration of Eq. (39) from t\ to ti results in 
Eq. (18) and thus both equations are equivalent. Therefore, the consistent time metrics 
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can be obtained, for example in the £ direction, by replacing • iy in the right-hand side 
of Eq. (34) with the time average, Vs</At, of Eq. (16) as, 



Note that £f defined by Eq. (40) contains all information about the movement of a cell 
surface, such as translation, rotation and deformation. In contrast, Eq. (34) is a sim- 
ple product of surface area and velocity of cell centroid and thus can represent only a 
translational motion. 

The free-stream subtraction technique will be useful for the rigid motion of the 
grid, because a rigorous evaluation of Eq. (40) is expensive computationally. Note that 
the subtraction is required only for the time-metric terms with the use of the free-stream 
capturing metrics in space. 5 


Concluding Remarks 

This paper summarizes the free-stream capturing techniques for the finite-volume 
(FV) and finite-difference (FD) formulations following a 1989 journal article by M. Vinokur. 
For an arbitrary motion of the grid, the FV analysis shows that volumes swept by all 
six surfaces of the cell have to be computed correctly. This means that the free-stream 
capturing time-metric terms should be calculated not only from a surface vector of a cell 
at a single time level, but also from a volume swept by the cell surface in space and time. 
The error introduced by conventional inconsistent time metrics is also shown. The FV 
analysis als o gives a guideline to construct free-stream capturing metrics in space and time 
for the FD formulation by regarding an FV cell as an FD mesh. The discretized geometric 
conservation law is shown to be a necessary condition but not a sufficient one. 
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Fig. 1 Geometry of a hexahedral cell. 
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Fig. 2 Volume swept by a surface. 
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